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04 ' We present a detailed semiclassical study of the propagation of a pair of optical fields in resonant 

^— , media with and without adiabatic approximation. In the case of near and on resonance excitation, we 

■ show detailed calculations, both analytically and numerically, on the extremely slowly propagating 

probe pulse and the subsequent regeneration of a pulse via a second pulse from the coupling laser. 
5— ( , Further discussions on the adiabatic approximation provide many subtle details of the processes, 

, ^ . including the effect on the pulse width of the regenerated optical field. Indeed, all features of the 

optical pulse regeneration and most of the intricate details of the process can be obtained with 
the present treatment without invoking a full field theoretical method. For very far off resonance 
excitation, we show that the analytical solution is nearly detuning independent, a surprising result 
that is vigorously tested and compared to numerical calculations with very good agreement. 
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Propagation of optical pulses in a resonant medium has long been an active field of research. The advances of 
. laser technologies have stimulated further development and understanding of the propagation of optical pulses under 
various conditions, in particular in the field of ultra-short and extremely high intensity optical pulse generation and 
propagation under adiabatic conditions. Early works in this field include the interpretation of the nonlinear wave 
mixing in terms of an adiabatic following model [1], counter-intuitive pulse sequence [2] and induced transparency 
[3]. Later, Eberly and co-workers [4] developed a theory for the propagation of optical field under the condition that 
atomic response can be treated adiabatically. Further development on the propagation properties of a pulse pair 
under the similar adiabatical conditions resulted in substantial understanding of the evolution of optical pulses in a 
resonant media. In particular, Grobe and co-workers first proposed the creation and recall of excitations distributed 
^ in media [5]. Recently, the field has experience a renewed interest and resurgence following the demonstration of 
extremely slow propagation of optical pulse both in ultra-cold and high temperature atomic gases [6] . Experimental 

■ studies on possible excitation of the atomic spin wave excitations with ultra-slow optical field, and the subsequent 
. ^ I "revival" of an optical field with frequencies near that of the ultra-slow field have promised future applications in 

■ information technology [7-9]. Theoretical works on such spin wave excitations and the subsequent "revival" of an 
Jh optical filed have also been formulated in a "dark-state polariton" theory [10], which basically is a field theoretical 

. . . reformulation of the adiabatic approximation of pulse pair propagation in a resonant medium originally formulated 
before [4,5]. These new treatments have provided some further understanding of the process, especially from the view 
point of atomic spin wave excitation and their relation to atomic coherence. 

Here, we present a study on the propagation and "retrieval" of ultra-slow optical fields. We investigate both 
analytically and numerically the propagation problem, with and without adiabatic approximation, for the cases of 
near resonance, on-resonance, and far detuned from resonance. In all cases, we show that the slow group velocity 
propagation and pulse "revival" is independent of the one-photon detuning by the probe laser, providing that the 
probe and coupling lasers are tuned to exact two-photon resonance and certain conditions are satisfied. That is, the 
results on off-resonance excitations are identical to the adiabatic theory formulation for the exact on-resonance case 
treated in [4]. The main contribution of our work, however, is the detailed analysis, both analytical and numerical, 
on the extremely slow propagation of the optical field, the corresponding atomic response, and the "revival" of an 
optical field when a second time-delayed coupling pulse is injected into the medium. These results, to the best of 
our knowledge, are not contained in the original studies by Eberly and co-workers and all the subsequent studies, 
including all recent theoretical works, on the subject. In addition, we show that all the experimental observables can 
be well predicted with the usual treatment combining the classical electrodynamics and the three-state model, without 
invoking a full field theoretical methodology. This is not surprising, since when a detailed handling of spontaneous 
emission is not critical and the fields are not extremely weak, a mean-field approximation to the quantum treatment of 
the electromagnetic field is valid. In Section II, we first present an adiabatic theory for the case where non-vanishing 
one-photon detuning exists. This treatment is very similar to the original treatment given in [4] except that the latter 
deals exclusively with the case of on resonance excitation. In Section III, we examine the propagation of a pulse 
pair and give two examples on how to obtain the evolution of the optical fields. Section IV constitutes the main 
contribution of the present study. We discuss in detail the propagation of the pulse pair in the adiabatic limit and 
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deduce analytically substantial insight and subtle understanding on the slow propagation of the probe field, the optical 
field-atomic spin cohenrence conversion, and the characteristics of the "retrieved" optical field. These results are not 
available in any of previous studies, and we believe they are very useful in providing a complete picture for the physics 
involved. In Section V, we compare predictions and estimates derived in Section IV with numerical simulations. 
Section VI is devoted to the case where the probe laser is tuned far from resonance, i.e. a Raman excitation scheme 
is employed. In this case, a different set of conditions are required for the adiabatic theory to be accurate. We show, 
however, that when certain conditions are satisfied, the solution is independent of the detuning, a surprising result 
that does not, to our best knowledge, exist in literature. Finally, the results are compared with numerical calculations 
for vigorous testing. In Section VII, we present a conclusion for our study. 



II. PRELIMINARY 



In this section, we examine the propagation of a pair of optical pulse in resonant media under the conditions where 
the adiabatic approximation can be made. We present here an expanded investigation based on the early treatment 
given in [4]. The difference is the inclusion of the non- vanishing one-photon detuning since the original treatment 
assumed exact on one-photon resonance excitations. We will show that this non- vanishing one-photon detuning leads 
to correct predictions about pulse propagation, coherent population transfer and storage, and the regeneration of an 
optical field when a time delayed second pulse is provided from the coupling laser. It should be pointed out that all 
theoretical studies published recently on this subject also assume exact on resonance tuning by the probe laser as 
originally assumed in Ref. [4] . Readers should consult a series pioneer studies on the on-resonance case given in [4] . 

Consider a A system as depicted in Figure 1. We assume that a probe laser {Ep, frequency LUp) is tuned on or near 
resonance with the |1 >^ |2 > transition, and a coupling laser {Ec, frequency ojc) is tuned so that exact two-photon 
resonance between states |1 > and |3 > is achieved. Assume atomic wave function of the form 

\^iz, tr) >= aie-''^i*|l > +a2e-''^^*\2 > +a3e-''^^*\3 >, (1) 

we thus find the following atomic equations of motion 

^ = i%e-"''''a2, (2a) 

^ = iQ.le^^-'ai+iQ.y^''as + i(5 + i'^^a2, (2b) 

^ = m,e-"''^a2, (2c) 

where (ij and 7j are the jth amplitude of the atomic wave function and decay rate, respectively. Qp and Qc ai'C! half- 
Rabi frequencies, given by ^* = D2iEp/{2h), ft* = D2zEc/ {2h) with being the dipole moment of the relevant 
transition. Also, tOp = u)2 — oji + 5 cjc = u!2 — cjs, and the quantity tr = t — z/c is the retarded time which is the 
combination of t and z that the laser field amplitudes depend on in vacuum. 

In order to correctly describe the propagation of the optical pulse, the atomic equations of motion (2) must 
be simultaneously solved self-consistently with Maxwell equations. In the limit of plane waves and slowly varying 
amplitudes, the positive frequency part of these fields satisfy 



dz 



= iKi2TAlA2, (3a) 



-^j = iKs2TA*A2, (3b) 

where the amplitudes Ai, A2, and A3 differ by position dependent phase factors from oi, 02, and 03, and we have 
introduced notations K12 = 2TrNu)p\Di2\'^ / (he) and K32 = 27rA^a;c| £'32 P/(?ic) with N being the concentration. 

In order to evaluate the polarization terms on the right hand sides of Eq.(3a,3b), it is necessary to obtain the 
amplitudes of the atomic wavehmction. The problem is intrinsically difficult because both Eq.(2) and (3) must 
be solved simultaneously. A great simplification can be achieved if one invokes the adiabatic approximation while 
evaluating the atomic wavefunction. The wavefunction obtained with the adiabatic limit will then be used to calculate 
the optical field in a self-consistent manner. This is the spirit of all adiabatic approximation based optical propagation 
problem and will be the method used in the present work. 
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In an adiabatic treatment wc introduce Ai = oi, A2 = e''^'p^a2, and ^3 = e^'^^^aa in order to remove the z dependent 
phase factors. The adiabatic eigenvalues of the equations of motion can then be obtained from the characteristic 
equation given by 



X'-{5 + i^) A2 - + |J7,|2)A = 0. 

They are 



Ao =0, (4a) 
A. = (^U,kp + W+(^V, (4b) 



As can be seen, the inclusion of the non- vanishing one-photon detuning changes the eigenvalues. For sufficiently small 

S, this condition reduces correctly to that given in [4] . 

The adiabatic condition requires that the eigenvalue XqT always differ from A±r by an amount that is large compared 
with unity. When \5t\ » \ilc\ and |i5r| >> \^p\, however, additional conditions are required. In fact, we require 
that the pulse length of the coupling laser be much longer than that of the probe laser, and |f7c'''P/|<^''"| >> 1 be 
satisfied. Therefore, during the entire pulse of the probe laser, the two-photon transition between |1 > and |3 > is 
shifted well outside the bandwidth of the probe laser. In another words, the overlap of the probe pulse band width 
with the resonance must be avoided. 

We now focus on the eigenvector corresponding to Aq = 0, since it is this adiabatic state that corresponds to the 
population all being in state |1 > before the laser pulses. In the lowest order of the adiabatic approximation Eqs.(2a) 
and (2c) imply that 02 = 0. Eq.(2b) then yields 

Ai = -^As. (5) 
Since, for a closed system, |aip -h |a2|^ + |(i3p = 1 we can write 

^\n,{z,trr + \np{z,tr)\-'' 

A2{z,tr) = 0, (6b) 
As{z,tr) = - , "^^^'^''^ (6c) 

It is important to realize that solutions (6a-6c) are not sufficiently accurate for many propagation problems. This is 
because that even though numerical solutions show that A^ is much smaller than the amplitudes of the other states, 
it cannot be exactly zero without assuming the laser pulses are propagating at exactly the vacuum speed of light. In 
many situations, however, this assumption of the speed of laser pulses is not valid. This is particularly true in the 
case of sufficiently high optical density and very near a strong resonance excitation where significant modification of 
the propagation velocity is anticipated. In order to accurately account the propagation effect in the regime where 
significant dispersion effect is encountered, one must seek the next order correction to the lowest order approximation 
where the small quantities such as A^ becomes critically important. To do this, one can use either Eq.(2a) or Eq.(2c) 
to derive a second approximation for A^-, thereby obtains the ffist order correction to the lowest order adiabatic 
approximation. Substitute Eq.(6a) into Eq.(2a) or Eq.(6c) into Eq.(2c) we obtain 

i dAx _ i d f fll \ 



and 
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Within the adiabatic approximation and the slowly varying amplitude approximation one can show that Qp(z,tr) 
and fic(^,tr) can be chosen to be real. It is then trivial to show that two forms for A2 as shown in Eqs.(7a,7b) are 
self consistent. It is interesting to note that even though having 6 non-zero makes the criteria for the validity of the 
adiabatic approximation different, it does not change the solution for the state amplitudes. 

As usual with the adiabatic approximation used in a resonance situation, we must have fig already strong when 
Sip starts to build up. This means that f2cT >> 1, where r is a measure of the time over which ilc, or flp change 
significantly. One scenario where the adiabatic approximation would hold through the whole laser pulse is if the two 
lasers peak at the same time, but the pulse length of the coupling laser is much longer. This feature of the pulse 
lengths and the requirement \0,ct\ » 1 through the pulse length of the probe laser will make results based on the 
adiabatic approximation quite accurate. In essence, one must keep both |A+t| >> 1 and |A_t| >> 1 throughout the 
pulse of the probe laser, e.g. there is a wide gap between Aq = and the other eigenvalues so that curve crossing is 
avoided throughout the entire pulse length. We will now investigate the use of the adiabatic approximation in the 
propagation problem. 



III. THE PROPAGATION PROBLEM 



In order to treat the propagation of the laser pulses through a resonant medium, the polarization of the medium 
as a function of position and time must be obtained. Using the atomic wave function obtained in Sec. II, we get 

P{t) = N < ^{z,tr)\D\^{z,t) > 

where N is the concentration. For the positive frequency parts of the polarizations at the probe and coupling laser 
frequencies, we have 

P+ = NDi2e'''"-'''^'''AlA2, (9a) 
P+ = A^£)32e''==-"^=*^^^2. (9b) 
Using the first order correction to A2, i.e. Eqs.(7a,7b), we find, for the amplitude of P+ 

1 d I SI* \ 

1 d f O* \ 

P^^ = ^^^Z2 . T7- , . (10b) 

Substitute these polarizations into Eq.(3a,3b), we obtain 

'dsi*p\ _ K12T d ( m; \ 



( dSl* \ K32T d I tSI* 



(lla) 
(lib) 



Equations (11) are the key equations describing propagation of a pulse pair within the adiabatic approximation. By 
using the adiabatic approximation for the atomic response one needs to solve only two nonlinear partial differential 
equations, instead of solving five nonlinear partial differential equations Eq.(2a-2c) and Eq.(3a,3b). 

To further illustrate the effectiveness of the adiabatic approximation, we now examine two special cases. 

We proceed by first noticing that the quantity 

i^ = L^ + L_£L^ (12) 
K12 K32 
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represents the sum of the photon fluxes at LUp and ujc divided by the concentration of the medium through which the 
waves propagate. Differentiate F with respect to z while holding tr fixed and make use of Eq.(ll), we obtain 



o , ' . " ■ ) 0. (13) 

OZ \ Ki2 K32 / 

Thus, F does not depends explicitly on z. This permits us to evaluate F by evaluating it at z = 0, the entrance to 
the atomic vapor cell, e.g. 

F(z,t)=F(z = 0,t)= I ' ^1 + I ' . 14 

K12 K32 

Thus, whenever \ilc{z,tr)\'^ / K32 + \^p{z, tr)\'^ / K12 occurs it can be replaced by F{tr), as determined in Eq.(14). An 
exact relation for this quantity can also be directly derived from Eq.(2a-2c) and Eq.(3a,3b). One obtains 

dz \ K12 K32 / dtr 

This relation indicates that anytime |^2p is very small and slowly varying, and the decay rate of this state is not too 
fast, this total photon flux is close to depending only on tr- The lack of dependency of this quantity on z when the 
full set of equations is solved numerically is an important test for the validity of the adiabatic approximation. 



Case 1. K12 = K32 

In this situation the quantity |r2p(z,tr)P + \i^c{z , U)]'^ = Ki2F{tr). We now use this fact in Eqs.(lla,b) to obtain 
an analytical approximation that can be compared with numerical solutions to Eqs.(2a-2c) and (11a, b). 

Since -y/|f2cP + |f^pP does not depend on z when tr is held flxed, dividing both sides of Eqs.(lla,b) by this quantity 
yields 



d I ^*pT \ _ K12T d I %T 



d ( ri*r \ K12T d ( VtlT 



dz y ^|n,T|2 + |r!pT|2 ) |r2cT|2 + |Qpr|2 dtr It \ ^|f]er|2 + |f2pT|2^ 

A few substitutions make an analytical solution to these equations relatively obvious. Let 



In terms of these dimensionless quantities Eqs.(15) become 



du dv 

where general "travelling wave" type solutions are immediately obtained as 



(15a) 
(15b) 



W„ = =, (16a) 

Wc = """^ (16b) 

V{tr) = f^'^ (|l^e(0, C)t|2 + \^p{Q, C)Tp) d {^-^ , (16c) 

^(-^) = / ni2Tdz . (16d) 
Jo 



dWc dWc „ , 
+ ^ = 0, (17b) 



Wp = Fp{v-u), (18a) 
Wc = F,{v-u). (18b) 
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This special case provides a very useful test case for examining the validity of the adiabatic approximation for the 
"trapped" light problem. The functions Fp and are easily determined by evaluating at z = (therefore, u = 0) and 
using the fact that the laser fields are known as a function of t and hence v. When the second coupling laser pulse 
is sent into the medium after a time delay, the predictions about the "revival" of the probe laser are all contained in 
this solution. The tabulation of Fp remains the same as long as only a second coupling laser is sent into the medium 
after some delay time. 



Case 2. |fij,(0,t)r| « \Clc{0,t)T\ with QpT « 1 

In this limit there can never be a significant population in state |3 >. Therefore, to a good approximation we have 

il.c{z,tr) = flc{0,tr). That is, the coupling laser is not effected very much by the resonant medium and it propagates 
at speed c without distortion. Thus, in this limit we only need to solve Eq.(lla) which becomes 



d ( n;T \ n^^T d ( n;T 



dz \ |J7c(0, tr)T\ J \nc{0, tr)T\^ dU / T \ \VLc{0, ir-)r| 

If we make the following dimensionless substitutions 



(19) 



v{tr) = f^'^ (|Oe(0,i;)r|2) d rA , (20b) 



u{z) = / Kl2Tdz , (20c) 



we then find that the solution for Wp is 

Wp = Fp{v - u), (21) 
with the function Fp being determined by using z = at the entrance of the cell, 

where v{t) is given in Eq.(20b), with = t at z = 0. A table of Fp as a function of its argument can be made, which 
will be used when z is no longer zero and the argument is ?;(tr) — u{z). Once we have tabulated Fp, the probe field 
at different z can be trivially calculated using the known coupling field, e.g. 

n;{z, tr) = |l^c(0, tr)\Fp{v{tr) - u{z)). (23) 

Again, the tabulation of Fp is determined completely at all later times by Eq (22) - at least for the case where only a 
time delayed coupling laser pulse enters the material at z = 0. 



IV. ADIABATIC APPROXIMATION APPLIED TO OPTICAL FIELD REGENERATION 



In the precious sections, we have provided detailed descriptions and examples on how adiabatic approximation 
works. In this section, we will discuss various aspects of the adiabatic approximation, especially the validity in the 
context of probe pulse "revival". These discussions provide many subtle understandings on the process which, in 
our view, are not discussed either in the recent works on "probe pulse revival" and the original studies of adiabatic 
approximation. This section, therefore, constitutes the main contributions of the present study. 

The validity of the adiabatic approximation depends on more than just the slow time dependence of the amplitudes 

of the probe and coupling lasers. For instance, the approximation will fail in dealing with the solution for the atomic 
response if the two pulses are sent simultaneously and the pulse widths and shapes are the same, even with both 
fields being turned on slowly. In this situation, at early times substantial populations will be promoted to state 
1 2 >, in contradiction to the adiabatic solution which, to the lowest order, predicts a zero population for this state. 
What is needed is to have the coupling laser built up and establish laser induced transparency before the probe laser 
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starts to become intense. This could be done by cither time delaying the probe laser relative to the coupling laser, 
or by making the coupling laser pulse length much longer than that of the probe laser and by having them peak in 
intensity at the same time. In the following treatment we will consider the latter situation. We therefore will choose 
the characteristics of ilp{0, t) and ^0(0, t) that would be conducive to making the adiabatic approximation work well, 
and demonstrate that this is indeed the case. We emphasize that our choice of the form of the laser fields is only one 
of many workable assumptions on the acceptable forms of the laser pulses that will make the adiabatic theory work 
well. 

Let us consider the functional forms for ilp{Q,t) and f7c(0,i) as 

Q.p{Q,t) = fipoe"^*/^)', (24a) 
Oc(0,t) = J7cO (e-0-2(*/^)' + ije-o-^f*/^-^")') . (24b) 

With this functional form, the Rabi frequency due to the coupling laser rises up to produce laser induced transparency 
before the probe laser becomes large enough to have any effect, an important consequence of having a significantly 
longer pulse length for the coupling laser. Without this characteristic of the coiipling laser, the probe laser would 
produce a population in state |2 > before the transparency was induced. The production of population in |2 > would 
mean that different adiabatic states have already been mixed and we have a failure of a strictly adiabatic picture 
in which only a single adiabatic state persists throughout the pulse. In Eq.(24), fipo and l^co are real constants 
characterising the peak amplitude of the two half- Rabi frequencies before the pulses enter the resonant medium. The 
parameter r is a measure of the pulse length of the probe laser, R is the ratio of the Rabi frequency at which the 
coupling laser recurs to its initial amplitude, xq = td/r is the value of tr/r at which the peak of the coupling laser 
recurs. 

We will begin by pointing out how the adiabatic approximation can be used to understand what happens for 
R = Q. That is, for the case where there is no recurring coupling laser pulse. We take fipoT = ^cot = 20, 
K12T = K32T = 700cm^^, and 72T = 0. We will soon see that in this numerical example the group velocity is 
sufficiently small so that the coupling laser dies away before the probe pulse can propagate through the vapor cell. 
When the coupling laser begins to die out, the intensity of the two lasers becomes proportional to each other. This 
is because when the coupling laser is beginning to become small, the following relation is appropriate (see Eq.(16c)) 

v{tr) |l^corpyW2+ |f^poTpyV2- (25) 

That is, almost all of the area under the two Gaussian pulse shapes has already been included in the integration. 
Also, |f2p(0,tr)|^ has long been very small compared with |f2c(0,tr)|^- Thus, for tr/r > 3.0 we have, as a very good 
approximation 

%{z,tr) ^ ^|f2c(0,ir)|2 + |Op(0,tOI^^;>((l^^cOT|^v/W2 + |Opor|2v/772) - M2T^). (26) 

(In the case of an inhomogeneous medium fi;i2r2; would be replaced by an integral over the z dependent K12T.) At 
such a late time during the pulse, since the argument of Fp depends only on the z coordinate, therefore the time 
dependence of il* is obviously exactly the same as that of y/\flc{0,tr)T\'^ + |Op(0, fr)Tp. When |flc(0,tr)| is several 
times larger than |rip(0,tr)|, as always is at such late times at z = 0, this means that fl*{z,tr) has the same time 
dependence as il*(0; tr) at siich late times. The fact that the two pulses have the same time dependence as they begin 
to become very weak is the reason that populations are left in each of the states |1 > and |3 >. If one of the laser 
fields died out much more rapidly than the other there would only be population left in only one state. 

^Prom Eq.(26) we see that if we choose z such that 

then the argument of Fp will be the same as at z = and tr = t = 0- With this argument, the value of Fp is 
Fp = fipo/i/l^coP + l^^pop- This is the largest value that Fp takes on. Thus, at this depth into the medium the value 
of 1^3(2:, tr)| matches its largest value at 2; = 0. This population persists at large tr until very slow collisional effects 
either destroy the coherence left behind in states |1 > and |3 >, or until collisions mix the populations of the two 
states. This long persistence of a coherent mix of populations in states |1 > and |3 > is what leads to the "revival" of 
an optical field with frequency very close to that of the original probe laser when a coupling laser is re-injected into 
the medium after some time delay [11]. The number of atoms left in |3 > is (within the adiabatic approximation) 
equal to the number of photons in the original probe pulse (recall that we have chosen parameters so that the probe 
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pulse docs not penetrate the entire medium). The creation of a coherent mixed stat is also what makes this situation 
interesting from the point of view of quantum computing. 

Next, we investigate the situation where i? > 0, i.e. a second couphng pulse enters the medium at a delayed time. 
Let flcOT = 20, flpoT = 5, 72T = 0, K12T = K32T = 200cm~^. From the above consideration, we expect lA^l to 
have its largest value at large U when z = 2.86 cm. In Figure 2 we show a contour plot of 1^3(2;, tr)| based on our 
adiabatic theory. The predicted IA3I is indicated by the line of constant color leading from tr = and z = out to 
the horizontal path at z ~ 2.86 cm, as expected. In Figure 3 we show a siufacc plot of the probe field Rabi frequency 
QpT as functions of tr/r at 2; = 3 cm. The long asymmetric tail of Qp is a consequence of the behavior described in 
Eq.(26). In the region between 2.5 < tr/r < 5 the ratio of the two half- Rabi frequencies is close to constant, averaging 
aroimd 0.24. This ratio is also close to the adiabatic approximation for ^3 since |rJpp << |r^cP- In a different set 
of parameters, we choose flcoT = ^poT = 20, 72r = 0, R = 2.923, and K12T = K32T = 700cm~^. In this example, our 
theory predicts that the peak value of \As{tr, z)\ before the second coupling laser pulse enters the medium will occur 
at z = 1.159 cm. Indeed, as can be seen clearly from Figure 4 that 1^3(2, t^)! takes on the value l/\/2 at a depth 
of about 1.1-1.2 cm into the medium, as predicted. Figure 5 depicts the corresponding probe field predicted with 
adiabatic theory for this set of parameters. 

We now illustrate how one can take a set of initial conditions and very simply determine the characteristics of the 

"revived" probe pulse at the point it emerges from the atomic vapor. 

First, note that with our choice of pulse characteristics and during the recurring coupling laser pulse 

viU) = l^corpyf (^l + er/(^y|i^^jj , (27a) 

S = |r!,oTpyW2 + |f^poT|2v/V2. (27b) 

The reader will recall that Fp was determined from the functional dependence of the probe and coupling laser pulses 
at 2; = 0. In particular, Wp{v{t)) = —A3{0,t), where ^3(0,*) is determined from Eq.(6c) and the pulse characteristics 
at 2; = 0. In addition, the following properties of the function Fp are very useful: 

Fp(0) = 0, (28a) 
Fp{S) = 0, (28b) 

Fp{S/2) = Fp{vm = "fj (28c) 

As will be seen later, we will invoke these properties when z = Zm is reached, where Zm is the value of z at the end of 
the medium, to estimate the time interval over which the "revived" probe pulse exits, as well as the peak amplitude 

of fl*{Zm,tr). 

Before proceeding further, we will point out some obvious properties that must be satisfied by the parameter, R, 
that appears in Eq.(24b). We first note that R must be large enough so that v{tr) — Ki2TZm > 0. That is, R must be 

large enough so that the group velocity becomes large enough to allow most of the regenerated photons to exit the 
cell before the laser induced transparency ends. This group velocity depends on position and time as (see Eq.(15)) 

^ = 1 + '^i^cr^ . . 

\np{z,tr)T\^ + \nc{z,tr)T\'^' ^ ' 



We first write down the condition for a time tri at which the "revived probe pulse" first reaches Zm- This is the 
earliest time at which the argument of Fp goes from being negative and becomes zero. At this point -Fp(O) = 0, but 
at later times it will become non-zero. This time is determined by 

V{tr\) - Ki2TZm = 0. (30) 

At time t^m the value of Fp{S/2) will be equal to —^3(0, 0). This is the maximum value Fp takes on with the set of 
parameters chosen. This time is determined by 

viUm) - Ki2TZm = S/2. (31) 

Finally, at the time tr2 at which the "revived probe pulse" completes its exit from the medium. At this time, we have 
Fp{S) = 0. Therefore, 

V{tr2) - Ki2TZm = S. (32) 
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Using Eqs.(27a,b) we can rewrite Eqs. (30-32) as 

:{trl - td) 



erf 



2(ki2T2to - S) 



erf v/2/5 



erf v^2/5 



T 

{tr2 — td) 



2{Ki2TZm - S/2) 
2{Ki2TZra) 



1, 
1. 

1, 



(33a) 
(33b) 
(33c) 



which can be further condensed into 



where 



erf 



/3 



{tr-td)\ __ a-fP 

i?2 



2Ki2TZm 

|ficor|2v/W2' 
25* 



- 1, 



with / = 1 to determine tri, f = 0.5 to determine trm, and / = to determine tr2- Correspondingly, we obtain the 
peak value of ^l*{zm, trm)'T as 



-(trm-td) /(5r )_ 



n 



p07 



(34) 



•pO^ 



As an example, let K12T = k^2T = 200cm ^, z„j = 8 cm, ilco''" = 20, flp^r = 5, i? = 4, and 72T = 0. With these 
parameters, we find a = 2.854598, a/i?^ = 0.1784124, /3 = 2.0559017, and P/R^ = 0.128494. Thus, {tri - td)/T = 
—2.19, {trm — td)/T = —1.767, and {tr2 — td)/T = —1.50. The value of ^p{zrmtrm) is thus estimated as 



%{zm,trm) = 4(20) exp(-0.2(1.767)2) -^_= 



10.39. 



The time trm corresponds to the largest value of Fp, but not quite to the maximum of fl*{zm, tr). The actual maximum 
occurs closer to -1.74 with the maximum value of 10.44. However, our simple estimate usually does very well if R is 
large enough to easily let all of the regenerated probe laser photons escape. 

We next note that if R is chosen to make the argument of Fp exactly the same as its value at z = and tr = t = 
at the time tr = td, we then have 



1 + 



Rf 



— Ki2TZm = 



Thus, 



|^^pOrpV^+ \^cOT\^V5^{l + R'^) = Ki2TZm. 



(35) 



(36) 



By choosing the value of R such that this equation is satisfied, one achieves a situation where the peak value of 
^p{zm, td) is given by (see Eq.(34)) 



QpOT 



(37) 



The charac;tc!ristics of this regenerated field can be estimated as follows. We first determine the value of tr such 
that the right hand side of Eq.(32) is zero or the full area under square of the incident probe and coupling lasers. 
This yields a range of time over which the probe rises from zero and returns to zero at the exit of the cell. One half 
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of this time interval is close to the full-width at half-maximum of the exiting probe laser beam. We then note that 
our condition on R yields 



Using these considerations we estimate the full-width at half-maximum pulse length, in the unit of the original probe 
pulse length r, to be 



_ 1 



ft 



cO 



(39) 



If the original probe laser has only penetrated a small fraction of the thickness of the medium when the coupling laser 
pulse has passed by, then R will turn out to be much larger than unity. This means that the width of the regenerated 
field is generally much smaller than the width of the initial probe pulse. Correspondingly, the bandwidth of the light 
will be larger. An immediate conclusion is thus available: the regenerated pulse is not the replica of the original 
probe pulse. It is coherent, to an extended degree, with the coupling laser because of the stimulated nature of the 
regeneration process. Therefore, the notion of "coherent storage" of the information carried in the original probe field 
is not accurate [7-9]. These features can be seen from Figures 3 and 5 which show f2* as a function oitr/r and z, as 
well as the time dependence of the regenerated field at the exit of the medium. In particular, as can be clearly seen 
from Figures 3 and 5, there is no optical field left in the region between the time when the original probe pulse and 
the companion coupling pulse enter the medium and the time when the second coupling pulse enters the medium. 
Therefore, photons arc neither "stored" nor "stopped" . All photons from the original probe pulse arc converted to 
the population of the state |3 >. These figures also show that in "reviving the probe pulse" every regenerated photon 
comes at the expense of flipping population from state |3 > to state |1 >. When the population of |3 > has been 
exhausted, there can be no further photon generated. Thus, in cases where the exiting regenerated field has a peak 
intensity much larger than the initial probe beam, the pulse width of this exiting pulse must be necessarily and 
correspondingly narrower. This is why the peak intensity of the regenerated field can go as R?, while the width of 
the pulse tends to go as l/R?. 

V. COMPARISON BETWEEN NUMERICS AND THE ADIABATIC THEORY 

We will now present some examples showing the degree of agreement between numerical calculations and the 
adiabatic theory developed in the previous sections. To obtain the numerical results, we solve numerically the coupled 
Maxwell and time dependent Schrodinger equations 



dtr 



iOp^2, (40a) 



^) = i%A^+iniA:, + i{5 + i^^A2, (40b) 
dA^ 



dtr 



iOc^2, (40c) 

(40d) 



dz 
dz 



iKi2AlA2, (41a) 
iK32AlA2. (41b) 



Here, Ai,A2, and A^ are the same as in Eqs.(2a-2c) except for the removal of the z dependent phase factors. First, 
recall the example that we have studied before where flcoT = 20, flpoT = 5, 72r = 0, K12T = K32T = 200cm~^. 
The discussions based on the adiabatic approximation predicts that a population of l^sl = l/\/l7 = 0.242536 at a 
depth of z = 2.86 for tr/r > 6.0, as can be seen from Fig. 2. As a comparison, we show, in Figure 6, a contour 
plot of 1^3(2,^^)1 for the same parameters given in Figure 2 but using numerical solution of Eqs. (40-41). The flat 
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region in the middle, a region where both of the probe and couphng lasers have been extinguished, indicates that 
during this period of time the population in the two-photon state remains nearly constant, just as we have concluded 
before based on the adiabatic theory that a population is left behind in states |3 > and |1 > which persists for a 
long time. Correspondingly, we show, in Figure 7, a surface plot of the probe field resulted from the direct numerical 
solution of Eqs. (40-41) using the same parameters as in Figure 3. We note that there is no probe laser field between 
tr/r = 2 and the time of the injection of the second coupling laser pulse, the same period where the population in the 
two-photon state remains nearly constant. This indicates, as our adiabatic theory has shown, that photons are neither 
"stopped" or "stored" in the medium as claimed in Ref. [7-9]. As the original probe photons being slowed down, 
they are absorbed to produce the two-photon excitation with the above noted population in the state |3 >. When 
all probe photons have been absorbed, there will be no further increase of the population in state 3 >. One may still 
argue that the probe photons are "stored" in the form of atomic coherence involving the two-photon excitation, and 
later when a second coupling pulse is injected into the medium, the probe pulse is "retrived" from this persist atomic 
coherence. This is, however, not correct. As we have shown the regenerated pulse has very different time-spectra 
characteristics from that of the original probe pulse. Therefore, the concept of "retriving the probe pulse" is incorrect. 
Indeed, the Fourier transform spectrum of all the regenerated fields have shown both different frequency content and 
noise characteristics that were not there in the original probe pulse. A simple argument based on our result (see 
Eqs. (37,39)) is sufficient to show that the regenerated pulse is very different from the original probe pulse, that is the 
regenerated pidse contains information in the second coupling pulse. Therefore, can not be the faithful replica of the 
original probe pulse alone. 

Further comparison between the adiabatic theory and the numerical solution are presented in Figures 8 and 9 for the 
second set of parameters used to create Figures 4 and 5 which were obtained with adiabatic theory. Finally, we show 
a case where non-vanishing one-photon detuning is assumed. Figure 10 is surface plot of the probe field obtained 
from the adiabatic theory with 5t = 10, whereas Figure 11 is obtained by directly solving Eqs. (40-41) using the same 
parameters. All these comparisons show very good agreement between our theory and numerical solutions. 

The numerical calculations were also tested by making use of the fact that with 72 r = 0, the sum of the squares of 
the three state amplitudes should be unity. In all numerical examples described in this work, this condition of unity 
was preserved through at least seven significant figures if 72T = 0. Also, in all cases when the number of tr/r and z 
grid points were doubles for all figures shown, the values of the state amplitudes and half-Rabi frequencies repeated 
through five significant figures at all points where the functions were not near a zero. 



In this Section, we examine the far off resonance excitation, a case has never been studied in literature with the 
same adiabatic treatment described in previous sections. Specifically, we investigate the condition for Aqt to stay 
well separated from \±t. We will show analytically that the formulism for far detuned excitation exhibits detuning 
independent feature when a specific condition is satisfied, a surprising result that to the best of our knowledge has 
never been reported in the literature. 

It is easy to see the condition for the eigenvalues to be well separated when \5\ » \flc\ and |^| >> In this 
limit 



which could bec;ome very close to the eigenvalue of the dark state lamdao = at sufficiently large detuning. Since 
the length of the pulse due to the coupling field is chosen to be much longer than that due to the probe laser, this 
eigenvalue can only be well separated from the zero eigenvalue at times when the probe laser first starts to build up 
if |A_t| >> 1 at such times. Physically, what is required is for the ac Stark shift in level |3 > to satisfy 



VI. CASE OF FAR FROM RESONANCE 






T >> 1. 



s 
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This means that as the two-photon transition between |1 > and |3 > is driven, the final state, i.e. state |3 >, must 
be shifted out of exact resonance by an amount that is much larger than the laser bandwidth. (We have assumed 
that the natural width of level |3 > is very small, being due to coUisional relaxation in a very cold vapor. If the laser 
bandwidth is less than the width of the state |3 >, the latter will replace in the above equation.) In Figure 12, 
we first show a surface plot of flp{z, t) based on our adiabatic theory for zero detuning. Figure 13 is a corresponding 
plot where the full numerical evaluations of Eqs. (40-41) are carried out. As expected, the solution based on adiabatic 
theory agrees well with the numerical solution. We now compare Figures 12 and 13 with Figures 14 and 15 which 
are the surface plots of ilp{z,t) with large one-photon detuning St = 120. Figure 14 is based on our adiabatic theory, 
therefore, should be compared with Figure 12, whereas Figure 15 is the result of full numerical calculation, thus 
should be compared with Figure 13. From these figures, we therefore conclude that when the appropriate adiabatic 
condition is met, the solution to the problem is insensitive to the detuning. Indeed, when the conditions given above 
are satisfied, the adiabatic approximation will be valid for all smaller S as well, a surprising results that has not been 
reported in literature and could play an important role in the future experiments. 

VII. CONCLUSION 

In conclusion, we have shown in detail how the adiabatic approximation can be used to understand the propagation 
of a pair of optical pulses in a highly dispersive medium. Starting from the adiabatic solution of a three-level system 
interacts with two laser fields, we showed analytically that in the case where the one-photon detuning is small (i.e. 
< 6 < \^c\), a coherent optical field with the frequency very close to that of the original probe can be regenerated 
when a coupling laser pulse in re-injected into the medium at a delayed time. We have provided detailed analysis on 
the conditions and characteristics of the regenerated field, including the estimate of the pulse width. Our analysis has 
shown that the regenerated pulse has a very different width therefore cannot be a perfect replica of the original probe 
pulse. Therefore, the notion of "coherent storage of the original probe pulse" is an incorrect description. As expected, 
the regenerated pulse is coherent, to an extended degree, with the coupling laser but not the original probe field. 
Subsequent numerical calculations have confirmed these analytical results. We have further developed the adiabatic 
theory to the cases where the one-photon detuning is larger enough so that |0c| < S and obtained a different type of 
adiabatic condition and solution. For time varying probe and coupling fields, the highly nonlinear wave equation for 
the recurring probe field is then investigated numerically. It is shown that even with this large one-photon detuning, 
as long as one chooses the pulse envelope and sequence properly so that there is a coherence left in the system, then 
there will be a coherent optical field generated when the coupling laser is turned on only. 

We thank Professor R. Grobe for providing reprints. 
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Figure Captions 

Figure 1. Energy level diagram showing relevant laser excitations. The decay rate of the state |3 > is assumed to be very 
small and hance neglected. 

Figure 2. Contour plot of A3{z,t) based on adiabatic theory. Parameters used: QpT = 5, = 20, 72T = 0, K12T = 
K32T = 200cm-\ i? = 4, td/r = 11. 

Figure 3. Surface plot of Qp{z,t) based on adiabatic theory. Parameters used: flpT = 5, CIcT = 20, 72T = 0, /ti2T = 

K32r = 200cTO-\ i? = 4, td/r = 11. 

Figure 4. Contour plot of ^3(2, t) based on adiabatic theory. Parameters used: QpT = flcT = 20, 72r = 0, K12T = K32T = 
700cm-\ i? = 3, td/r = 11. 

Figure 5. Surface plot of ilp{z,t) based on adiabatic theory. Parameters used: flpT = ^lcT — 20, 72r — 0, K12T — K32T = 
700cm-S R = 3, td/r = 11. 

Figure 6. Contour plot of ^3(0,^) based on numerical solution of Eqs. (40-41). Parameters used: VlpT = 5, Q.cT = 20, 
72T = 0, K12T = Ks2T — 200cm~^, R = A, td/r = 11. This figure should be compared with the Figure 2. 

Figure 7. Surface plot of Vtp{z,t) based on numerical solution of Eqs. (40-41). Parameters used: VipT = 5, flcT = 20, 
72T = 0, K12T = Hi2T — 200cm^^, i? = 4, td/r = 11. This figure should be compared to Figure 3. 

Figure 8. Contour plot of Az{z, t) based on numerical solution of Eqs. (40-41). Parameters used: ilpT = ficT = 20, 72T = 0, 
K12T = K32T = 700cm~^, R = 3, td/r = 11. This figure should be compared with Figure 4. 

Figure 9. Surface plot of 0,p{z, t) based on numerical solution of Eqs. (40-41). Parameters used: flpT = = 20, 72T = 0, 
K12T = K32T = 700cm~^, R = 3, td/r = 11. This figure should be compared with Figure 5. 

Figure 10. Surface plot of ^lp{z,t) based on adiabatic theory with non-vanishing one-photon detuning. Parameters used: 

QpT = 10, = 20, 72T = 0, (5t 10, K12T = K32T 200cm-i, i? = 4, td/r = 11. 

Figure 11. Surface plot of Qp{z,t) based on numerical solution of Eqs. (40-41) with non-vanishing one-photon detuning. 
Parameters used: flpT = 10, = 20, 72r = 0, St = 10, K12T = K32r = 200cm~^, ii = 4, td/r = 11. 

Figure 12. Surface plot of flp{z,t) based on adiabatic theory with zero detuning. Parameters used: flpT = 10, CIcT = 40, 
72r = 0, 6t = 0, K12T = K;32r = 200cm~^, R = 4, td/r = 10. 

Figure 13. Surface plot of ^lp{z,t) based on numerical evaluation of Eq.s(40-41) with zero detuning. Parameters used: 
ilpT = 10, ^IcT = 40, 72r = 0, 5t = 0, k\2T = K32''" = 200cm~^, R = A, td/r = 10. This figure should be 
compared with Figure 12. 

Figure 14. Surface plot of ^p{z, t) based on adiabatic theory with large one-photon detuning. Parameters used: fipT = 10, 
VLcT = 40, 72T = 0, i5t = 1200, K12T = K32T = 200cm~^, R = A, td/r = 10. This figure should be compared 

with Figure 12. 

Figure 15. Surface plot of f2p(z, t) based on numerical evaluation of Eq.s(40-41) with large one-photon detuning. Parameters 
used: flpT = 10, CIcT = 40, 72T = 0, 6t = 120, K12T = K32T = 200cm~^, R = 4, td/r = 10. This figure should 
be compared with Figure 13. 
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